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A 20 Year Radio Light Curve for the Young Supernova Remnant 
Gl.9+0.3 
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ABSTRACT 

The radio source Gl .9+0.3 has recently been identified as the youngest known Galactic super- 
nova remnant, with a putative age of ^ 100 years. We present a radio light curve for Gl .9+0.3 
based on 25 epochs of observation with the Molonglo Observatory Synthesis Telescope, span- 
ning 20 years from 1988 to 2007. These observations are all at the same frequency (843 MHz) 
and comparable resolutions (43" x 91" or 43" x 95") and cover one fifth of the estimated 
lifetime of the supernova remnant. We find that the flux density has increased at a rate of 
over the last two decades, suggesting that Gl.9+0.3 is undergoing a 



1.22± °-Jg percent yr" 



period of magnetic field amplification. 
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1 INTRODUCTION 

The well known deficit of young supernova remnants (SNRs) in 
our Ga laxy has motivated many searches f or these objects (for ex - 
ample; iGreen & Gulilll984 lGravlll994al : iMisanovic et"ani2002h . 
Predictions based on extragalactic supernova (SN) rates suggest 
that there sho uld be around 40 SNRs younger than 2000 years 
in our Galaxy (Cappellaro 2003), and yet less than 10 have been 
identified tGree n 2004). The recent ide ntification of Gl.9+0.3 as a 
, young SNR with a ge ~ 100 years by [Reynolds et al.l ( l2008h and 
iGreen et"al] ( |2008|) . adds to this set, and makes it potentially the 
youngest known Galactic supernova remnant. 

Most Galactic SNRs detected at radio wavelengths are de- 
tected well after the initial sup ernova event (> 1 0"^ — 10^ years). 
Cassiopeia A (SN 1681±19: Ipesen et al.|[2006h is the youngest 
SNR of known age, as ~ 330 years. The early stages of radio su- 
pernova (RSN) evolution have been studied in a number of bright 
extragalactic sources; detailed studies include, SN 1979C in Ml 00 
dWeiler et alJ[T986) . SN 1980K in NGC 6946 jWeiler et alJll986l■ 
ll992l) ■ SN 1993J in NGC 3031 ( IWeiler et al.ll2007h and SN 1987A 
in the Large Magellanic Cloud ' Ball et alj|200ll ; Manchester et al.1 
l2002h. Ho wever, the oldest RSN of known age is SN 1923A 
teck et al]|l998 ). which leaves a critical gap in our knowledge of 
supernova evolution at intermediate ages of ~ 50 — 300 years. 
Hence the age estimate of ~ 100 years for Gl.9+0.3 makes it a 
useful probe of this period in SNR evolution. 

Gl.9+0.3 was first identifi ed as a potential you ng SNR based 
on its small angular size of 1.2'(|G reen &Gullll984 . Recent Chan- 
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dra observations bylR evnolds et al.l (l2008l) showed that the remnant 
had expanded by 16 ± 3 per cent between earlier VLA observations 
in 1985 and the Chandra observatio ns in 2007. This w as confirmed 
by VLA follow-up observations bv lGreen et al.l ( l2008h which gave 
an expansion of 15 per cent over 23 years (or an expansion rate of 
0.65 per cent yr^^). This led to an a ge estimate of ^ 100 years 
(or at most 150 years) for Gl.9+0.3. iGreen et al.l ( l2008h also use 
archival data from a variety of instruments (over frequencies 332 
to 5000 MHz) to suggest that Gl.9+0.3 has brightened by a rate of 
~ 2 per cent per year over the same period. 

We present 20 years of radio observati ons from the Molonglo 
Obse rvatory Synthesis Telescope (MOST; lMillslll98ll ; iRobertsonI 
fl99lh . These observations were carried out at constant frequency 
(843 MHz) and comparable resolutions (43" x 91" or 43" x 95"). 
They show that Gl.9+0.3 has been increasing in brightness by a 
rate of 1.22 ± [j'^g per cent yr~^ over this period. Our new es- 
timate has the advantage that th e measurements we re taken with 
the same instrument, whereas the lCreen et alj ( I2OO8I) estimate was 
based on observations from a range of instruments, compiled from 
the literature. 



2 OBSERVATIONS AND DATA 

The observations presented have been carried out over the last 20 
years with the Molonglo Observatory Synthesis Telescope. The 
MOST operates at a central frequency of 843 MHz and has a restor- 
ing beam of 43" x 43" esc \S\ where S is the c entral declination o f 
the field. Further technical specifications are in iBock etal.l ( fT99^ . 
Gl.9+0.3 is detected in a total of 30 archival fields covering 
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the period 1988 to 2007. These fields were observed as part of on- 
going projec ts; the Molonglo Galactc Centre Survey (1985-1991; 



Grav '1994b'), the Molonglo Galactic Plane Survey (1983-1994; 
Green et al. 1999 ), the second epoch M olonglo Galactic Plane Sur- 
vey (1996-2007; iMurphv et aLllOOTl) , and a Galactic centre mon- 
itoring program (2004—2007). Discarding five fields in which the 
data quality was poor, we have 25 observations of Gl. 9-1-0. 3. The 
field of view for the early observations was 70', which increased to 
163' after the MOST was upgraded in 1995. The restoring beam is 
43" X 95" for the pre- 1995 fields and 43" x 91" for the post- 1995 
fields. Therefore the frequency and resolution were comparable for 
all of our observations, making this a self-consistent dataset cover- 
ing approximately one fifth of the life of Gl. 9-1-0.3. 

The MOST data are reduced using a custom process described 
by iBock et al.l ( Il999h . The phase and flux density calibration are 
done by observing a set of calibrator sources before and after 
each 12 hour observation. These cal ibrators are strong, unresolved 
sources taken from the list given in Campbell- Wilson & HunsteadI 
1994j). Analysis of the calibrator sources by Gaensler & HunsteadI 
2000f ) showed the relative flux scale (as measured by the aver- 



age scatter in repeat observations of the non-variable calibrators) 
to be accurate to ar ound 3 per cent. A similar analysis done by 
iMurphv et al] ( |2007|) , found the scatter for non-calibrator sources 
observed by MOST to be around 5 per cent. 

The noise in individual MOST images is a combination of 
thermal noise and source confusion from both the main beam and 
sidelobes of the MOST. The noise level varies with Declination. 
Since the telescope is calibrated only at the beginning and end of 
each 12 hour run, there may be changes in the precise operational 
condition of the telescope electronics and feeds that have a small 
effect on the sensitivity of a given observation. The rms noise at 
the centre of a typical MOST image is ~ 1 mjy beam^^. However, 
most of the fields containing G 1.9-1-0. 3 have strong artefacts caused 
by the Galactic centre which is also in the field. Hence typical rms 
noise levels near the source are ~ 50 mJy beam^^ . 



3 ANALYSIS 

The MOST observations are not at high enough resolution to detect 
significant changes in the size or mophology of G 1 .9-1-0.3 with time. 
Gl. 9-1-0.3 is only slightly extended (~ 2 beams E-W) in the MOST 
images, as shown in Fig.[T]which shows a selection of our observed 
epoch s. Fig. [2] shows a M OST epoch (from 2007/06/19) overlaid 
on the lGreenetall ( l2008h VLA 4.8 GHz image (from 2008/03/12). 
When the VLA image is smoothed to the MOST resolution, there 
is good agreement in the observed morphologies. The rms noise 
in the 4.8 GHz image is ~ 0.012 mJy beam"\ and in the MOST 
image is ~ 30 mJy beam~^. 

The integrated flux density of the source was measured by first 
fitting a two dimensional second order polynomial surface to the 
surrounding region and subtracting it, and then fitting a double el- 
liptical Gaussian to the source itself. Fitting the background was 
necessary due to the large flux gradient in the image caused by the 
artefacts from the Galactic center. In most of the fields the source 
is sitting in a significant negative trough that could have negative 
flux values as large as —100 mJy. Although the VLA observations 
show that morphology of the source is clearly a shell, at the MOST 
resolution it can be modelled reasonably well with a double ellipti- 
cal Gaussian. This is demonstrated in the lower two panels of Fig.|2] 
which show a double Gaussian model, and the residuals after sub- 
tracting the model from the VLA image smoothed to the MOST 
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Figure 3. 843 MHz radio light curve for Gl. 9+0.3 from 1988 to 2007. Flux 
densities were calculated using a double Gaussian fit after polynomial sub- 
traction of the background (see text for details of the flux density and eiTor 
calculations). The solid line shows a least squares fit with gradient 0.015 
Jy yr-i and a flux density of 1.23 Jy on 2005/01/01. 



resolution. In the MOST images, the residuals after fitting were 
around the ~ 5 per cent level, which we have included as a com- 
ponent of our overall error. The resulting light curve is shown in 
Fig. [3] Errors on each point are calculated by adding the estimated 
error in the flux density scale (~ 5 per cent), rms noise (calculated 
in region local to the source in each map, typically 20 — 50 mJy 
beam^^ or ~ 5 per cent) and estimated fitting errors (~ 5 per cent) 
in quadrature. 

Fitting a power law (/ oc t™) to the light curve resulted in an 
exponent of m = 1.2 with large errors. A least squares fit of the 
light curve with a straight line gave a gradient of 0.015 Jy yr~^, 
with a reduced of 2.0 for 25 degrees of freedom. Since this 
does not allow an accurate estimate of our er rors, we have per- 
formed a bootstrap analysis, as described by lEfronl ( Il982l) . For 
our analysis we selected 25 000 random samples with replacement 
('bootstrap samples') of our flux measurements and fit each sample 
to form a probability distribution of the resulting parameters. From 
this we measured a best-fitting gradient and 68 per cent-confidence 
range of 0.015± Jj ^pj Jy yr~^. This analysis allows us to rule out a 
flat or decreasing flux density (a line with gradient ^ 0) with 99.95 
per cent confidence. We also constructed light curves for several 
field sources common to most of the epochs in our sample. There 
was no systematic change in flux density over this period of our 
observations. Comparison with our fitted flux value of 1.23 Jy on 
2005/01/01 gives a percentage increase of 1.22 per cent yr~^ with 
68 per cent-confidence range of 1.06 — 1.46 per cent yr"'^. 

Taking the age estimate of ~ 100 years from iGreen et"al] 
(2008), we have plotted the light curve of G 1.9-1-0.3 in Fig. |4] 
with a selection of other RSNe (SN 1923A, SN 1950B, SN 1957D, 
1979C, 1980K) and SNRs (Cas A, Kepler - SN 1504, Tycho - SN 
1572) for comparison. To convert flux density to luminosity, we as- 
su med a distance of 8.5 kpc for Gl. 9-1-0.3, following the argument 
of lRevnolds et alj ( l2008h that the high absorption towards Gl.9-1-0.3 
makes it unlikely that it is significantly closer than the Galactic 
Center. For plotting clarity, the points have been binned into three 
groups of measurements (~ 1990, ~ 1998, ~ 2005). 
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Figure 1. MOST images for a selection of our observations, before the polynomial background subtraction. Contours are at 50, 100, 200, 500 mJy beam ^ 
and the greyscale range is —200 to —350 mJy beam~^. The MOST beam FWHM is shown in the lower right of each image. 



4 DISCUSSION 

Only a fraction of SNe have detectabl e radio emission . Observa- 
tions of these RSNe (see, for example; lEck et al1l2002h show that 
radio emission caused by the interaction between the shock and the 
circumstellar medium of the progenitor star is first detected within 
days to months of the initial explosion. As the shock travels into 
regions of decreasing opacity, the radio emission brightens, peak- 
ing between days to years after the explosion. After reaching their 
maximum brightn ess, RSNe then show a power law decrease in flux 
density with time jWeiler et alj|2002h . This is demonstrated by the 
hght curves of SN 1923A, SN 1950B, and others, shown in Fig.|4] 
The radio behaviour o f RSNe is explained well by the models of 
IChevalieilfT983.ll998h . 

There are even fewer Galactic SNRs of known age. Cas- 
siopeia A, the youngest SNR of known age at ~ 330 years, Kepler 
(SN 1604) and Tycho (SN 1572) are shown in Fig.H The radio 
emission we detect from SNRs is due to the interaction between 
the shock and the interstellar medium. Hence the radio emission 
from a young SNR will increase once it has swept up enough of the 
surrounding interstellar mediu m. The timescale of this brightenin g 
is predicted to be ~ 100 years ( lGul]|l973l ; ICowsik & Sarkaill984h . 

Between the youngest SNR of known age (Cas A) and the 
oldest RSNe of known age (SN 1923 A), there is a gap in our ob- 
servational evidence which makes it hard to probe the period af- 
ter the fading of t he radio emission from the RSN and before the 
SNR switches on. lEck et alj j2002h conducted a search for radio 
emission from SNe of a range of known ages, but did not de- 



tect any of the SNe that occurred prior to SN 1923 (for example, 
SN 1885A and SN 1909 A). If the estimated age of 100 - 120 years 
for G 1.9-1-0. 3 is correct, then our measurements help constrain the 
flux density evolution in this intermediate time range. 

To ex plain the observed expansion rate measured by 
iRevnoldre t al. (2008) and the new light curve presented here, we 
need to reconsider some of the standard assumptions made when 
predicting the luminosity evolution of SNRs. The radio luminosity 
of a synchrotron source at a given frequency f is a function of 
the energy spectrum of the ultrarelativistic electrons, the magnetic 
field B present in the s ource, and the volume of the source V, as 
given m iLongaiil fT994h 

= A{a)VKB^+°'u-°' , 

where a is the spectral index, A{a) is a constant, and k is defined 
in terms of the electron energy spectrum per unit volume 

N{E)dE = KE-^^°'+'-'>dE . 

Hence for a given frequency 

If we make the assumptions that the expansion of the remnant is 
adiabatic, and that magnetic flux freezing is applicable, then the 
magn etic field strength decreases as B r^^ . Following iLongaiJ 
Il994i) we can derive that Vn oc r and so 

Now the time dependence of the luminosity is proportional to the 
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Figure 2. A comparison of one MOST epoch (2007/06/19) with a VLA 4.8 GHz image from lGreen eTaP <2008h . Top left: MOST contours (200, 400 and 600 
mJy beam~^) overlaid on VLA greyscale (range to 13 mjy beam~^). The MOST and VLA beams are shown to the lower right. Top right: The VLA image 
smoothed to MOST resolution is shown in contours (100, 150 and 200 mJy beam^^) overlaid on the original image. Bottom left: A double Gaussian model 
fit is shown in greyscale (range to 230 mJy beam^^), with the smoothed VLA contours overlaid. Bottom right: The greyscale shows the residuals after 
subtracting the double Gaussian model fit from the smoothed VLA image (range —3 to 3 mJy beam~^) with the smoothed VLA contours for comparison. 



time dependence of the source radius, since the spectral index 
a should not change during adiabatic expansion. Expressing the 
change in radius with time as a power law with expansion parame- 
ter m, r oc t"^ , then 

During an ideal free expansion phase r oc t, and in the Sedov 
phase r oc t° Usin g the spectral index of q = 0.63 given by 
[Reynolds et al.l ( |2008|) we would expect the luminosity to change 
as Ljy oc t~*'^ in free expansion, or cx t~^'^ in Sedov expan- 
sion. 

If we fit the measured expansion rate given by [Reynolds et al.l 

( I2OO8I) with a power law (assuming an age of 100 years in 2008), 
we get r oc t"'^^, implying L„ oc t~^'"' which is intermediate 
between free expansion and the Sedov phase. Using our new light 
curve data we get oc t. Clearly, to produce the observed increase 
in luminosity, either the magnetic field strength or the energy den- 
sity of the relativistic electrons (or both) must be increasing with 
time, rather than d ecreasing as assum ed above. A similar conclu- 
sion was drawn bv lOreen et all feOOSi) . This scenario is supported 



by what we know about older SNRs — extrapolating backwards 
from the observed magnetic field strength in Cas A, the assump- 
tion of magnetic flux freezing would resu lt in an implau sible mag- 
netic field strength in the progenitor star ('Longair"l994'). Hence at 
some period in its evolution, the magnetic field strength must have 
underg one a process of amplification. Simulations bv lJun & Jones! 
Il999j show that a possible mechanism for this amplification is 
Rayleigh-Taylor instability in the interaction region between the 
SNR shock and a surrounding cloud of interstellar gas. This hap- 
pens on timescales of ~ 100 — 300 years after the inital explosion. 
Our observations suggest that Gl. 9+0.3 has entered such a phase, 
in which kB is growing. 



5 SUMMARY 

Twenty years of observations with the Molonglo Observatory Syn- 
thesis Telescope show that the young supernova remnant G 1.9+0. 3 
has increased in brightness by 1.22 ± [[ ^g per cent yr~^ between 
1998 and 2007. This supports the estimate of ~ 2 per cent yr~^ 
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Figure 4. Comparison of the ~ 1.4 GHz light curve of Gl.9+0.3 with 
those of other RSNe and SNRs of known ages. The MOST flux densi- 
ties of Gl.9+0.3 have been scaled to 1.4 GHz assuming a spectral in- 
dex of a = 0.63. For clarity, the MOST data have been binned into 
three bins corresponding to the three groups of measurements (~ 1990, 
~ 1 998, ~ 2005). T he RSNe are SN 1923A, SN 1950B, SN 1957D (all 
from IStockdale et alj|2 006). SN 1979C tWeiler et al. 1991; M ontes et alj 
I2OOOI). SN 1980K iWe iler et al. 1986; Montes et al.ll 19981) and SN 1987A 
iStavelev-Smith et al.l l2007). The SNRs are Cas A (with the 1 965 - 1999 
fading shown, extended for clarity; Reichart & Step hensllOOOl) . Kepler (SN 
1604) and Tycho (SN 1572) (both from lGreeij|2004l) . 



made bv lGreenet This result, for what is potentially the 

youngest SNR of known age at ~ 100 years, confirms predictions 
that although most RSNe are decreasing in brightness, they will 
have to go through a period of brightening later in their evolution, 
to match the observed brightness of young SNRs in our Galaxy. 

As the next generation of radio telescopes revolutionise the 
study of RSNe and SNRs, archival data from telescopes such as 
MOST will have an important role to play, giving us insight into 
the evolution of these objects over long timescales. 
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